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Abstract — This paper presents modification of the conventional 
model reference adaptive control (MRAC) architecture in order 
to achieve guaranteed transient performance both in the output 
and input signals of an uncertain system. The proposed modifi- 
cation is based on the tracking error feedback to the reference 
model. It is shown that approach guarantees tracking of a given 
command and the ideal control signal (one that would be designed 
if the system were known) not only asymptotically but also 
in transient by a proper selection of the error feedback gain. 
The method prevents generation of high frequency oscillations 
that are unavoidable in conventional MRAC systems for large 
adaptation rates. The provided design guideline makes it possible 
to track a reference command of any magnitude form any 
initial position without re-tuning. The benefits of the method 
are demonstrated in simulations. 

I. Introduction 

The asymptotic behavior of adaptive systems is well re- 
searched during last couple of decades, and it is well known 
that the asymptotic tracking can be always achieved using 
so called Lyapunov redesign method. However, the transient 
behavior of the input and output signals can be very oscillatory 
with big excursions [13]. There has been great deal of efforts to 
modify the control architecture and the adaptive laws from the 
perspective of improving the transient behavior of the tracking 
error. The majority of these efforts led to high gain linear 
nonadaptive feedback [3], [4], [12], switching control law [8], 
[9] or to a parameter dependent persistent excitation condition 
[1]. 

In general, increasing the adaptation rate while reducing 
the tracking error magnitude (see for example [6]), generates 
high frequency oscillations and big overshoot in the control 
signal leading to possible actuator failures or excitation of 
unmodeled dynamics, which in turn can drive the overall 
system to instability. This shortcoming is common for the 
majority of existing adaptive control methods. Recently some 
results appear in the control community explicitly addressing 
the input signal transient behavior, which is very important 
from the point of view of the performance specifications of 
closed loop systems. 
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In [2] the authors proposed a new adaptive control architec- 
ture, called Li adaptive control, that can achieve close tracking 
of a given reference command both in transient and steady 
state by increasing the adaptation rate. The low pass filter 
introduced in the control channel prevents the high frequency 
oscillations in the control signal that closely follows the ideal 
one. However, the approach loses the possibility to have a real 
reference model to track, for which the control metrics can be 
specified. 

An alternative approach is proposed in [11] for a class of 
multi-input multi-output uncertain nonlinear systems to track 
a given reference model. Along with the tracking error it 
uses also its integral to guarantee the transient performance of 
both input and output signals. The control algorithm internally 
generates a low pass filter, thus preventing high frequency 
oscillations for the large adaptation rates. 

In this paper we take somewhat different approach and in- 
stead of modifying the control architecture or the adaptive laws 
we modify the reference model by feeding back the tracking 
error signal. Similar approach was used in [7] to improve 
robustness of the system with respect to disturbances which are 
bounded by the tracking error. Unfortunately the authors did 
not further investigate the properties of the control algorithm. 
We use the error feedback method from different perspectives. 
Our approach, called modified reference model MRAC or M- 
MRAC in short, is motivated by the fact that the initial large 
error in the control gains generates large transient excursions 
both in system’s control and output signals. Moreover, the out- 
put excursion in its turn generates oscillations in control signal. 
So the idea is to drive the reference model toward the system 
proportional to the tracking error, thus preventing the system’s 
attempt to aggressively maneuver toward the reference model. 
As the tracking error is vanished the modified reference model 
retains its original form. Therefore, the system asymptotically 
tracks not only the modified reference model, but also the 
original one. Moreover, the error feedback term determines 
the dumping in the control signal, increasing of which makes 
it possible to increase the learning rate for better transient 
performance without generating oscillations in the system. 
Meanwhile, very large error feedback gain may increase the 
overshoot with respect to original reference model. Therefore, 
both parameters need to be increased simultaneously for better 
transient performance, which can be made arbitrary close to 



the ideal performance of the known system. A design guideline 
is provided for the selection of the feedback gain relative to 
the adaptation rate. The proposed adaptive control method has 
uniform performance in the sense that the same tracking is 
achieved for different reference commands and from different 
initial conditions without re-tuning. 

The rest of the paper is organized as follows. Section II 
presents results from linear systems theory. Section III present 
the control design and defines the error signals. We analyze 
the asymptotic properties of the M-MRAC design in Section 
IV and the transient properties in Section V, where we also 
discuss the design specifics. A simulation example is presented 
in Section VI, and some concluding comments are given in 
Section VII. 

II. Preliminaries from Linear Systems Theory 

Consider a linear system y(s) = G(s)u(s), where G(s) 
is a proper stable transfer function. If u(t) G then the 
following inequality holds 

lli/WIUco ^IlflWIUJKOIkoo, (1) 

where g(t) is the impulse response of G(s) (see for example 
[14] (p. 110)). If u{t) G C 2 then we have [14] (p. 108) 

\\ymc 2 = \\G(s)\\ H J\u(t)\\c 2 , (2) 


On the other hand, it is possible to minimize a com- 
bination of ||Gi|| Woo and ||G , 2||7t 00 that results in u> inde- 
pendent minimizing (. One such combination is h((, n) = 
IIGiH^ + 1| G 2 1| , where n is a positive integer. Since h(C,) 
is homogeneous in u>, the minimizing value of ( is independent 
of u>. The minimum problem for « = 1 results in suboptimal 
values, which are close to the individual optimal values 

Co = 0.634, ||G 1 (s)||-ft = 4, IIGaWllwoo = - • (7) 

LU Z U! 

where c\ = 1.0198 and C2 = 1.4162. 

III. Reference Model Modification 
Consider a first order linear system 

x(t) = ax(t) + bu(t), x(0) = xo , (8) 

where x(t) G K and u(t) G 1 are the output and input of 
the system, and a and b are unknown constant parameters 
with known sign of b. To simplify notations and without loss 
of generality we assume that b > 0. In conventional model 
reference adaptive control (MRAC) framework, the objective 
is to design a control signal u(t) such that the output of the 
system tracks the output x°(t) of the reference model 

= -a m x 0 (t) + b m r(t), s o (0)=a;o, (9) 


Also, ||C?(s) H-Koo and ||g(f)||£ 1 satisfy the relationship (see 
Theorem 4.5 in [14]) 


ITOIIhoo < 


Cl < 2||G(s)|| Wo 


(3) 


Next, for each u>, we minimize the norms of the second 
order transfer functions 

1 


G l( S ) = -^ 2 

s z + 2Qujs + uo z 

r / \ _ s + 2 C<^ 

2( } s 2 + 2 (ojs + u 2 

with respect to (. Computation of the ||Gi(s)||-h 0 


|Gi| 


■Hex, — 


if C > 75 




if 


^ - C2 


(4) 

results in 

(5) 


minimum of which is reached at C = -75 and is equal to 


IGrH^ = The H a 

K 


V2 

norm G 2 (s) is computed to be 


IG 


2||Moo = 


^2^( v / l+87-2C 2 -l) ! 


c > Tip 
r < yi+c2 

8 — 2 


(6) 


minimum of which is reached at C = and is equal to 
II G2 ||hoo = notice that the functions || Gi || -^^ and 

II G2 ||?too reach the minimum with respect to ( at different 
points, that is there is no optimal value of ( independent of u 
that simultaneously minimizes both IIGiH?^ and || G2 - 
Instead, the whole interval < C < 75 represents subopti- 


mal values for £, which results in 

1 „ 1.0328 1.4142 

— < HGrll^ < 


<I|G: 


2IIK00 < 


1.4553 


where a m > 0 and b rn are chosen to meet the performance 
specifications, and r(t) is a smooth and bounded external 
command. We notice that this reference model results from 
the application of the ideal control signal 

U°(t) = k\x°(t) + fc2 r (f) , (10) 

where the ideal control gains k\ and fcjj satisfy the matching 
conditions 


bk\ = —a — a m , bkl = b m . (11) 

The reference model (9) always can be specified from the 
performance perspective, but the ideal control signal (10) can 
not be implemented since the parameters a and b are unknown. 
Therefore, in MRAC framework the adaptive version of the 
ideal control signal is implemented, that is 

u(t) = ki(t)x(t) + k 2 {t)r{t) , (12) 

where ki(t) and k 2 (t) are the estimates of the ideal control 
gains k\ and k 2 , and are generated by the adaptive laws 

h (t) = — 7x(f)e°(f), ki(0)=k 10 

k 2 {t) = -7 r(f)e°(f), fc 2 (0) = fc 2 o , (13) 

where 7 > 0 is the adaptation rate, and e°(t) = x(t) — x° (t) is 
the tracking error. Introducing the estimation errors as k\(t) = 
ki(t) — k'i and k 2 (t) = k 2 (t) — k 2 , and taking into account 
the matching conditions (11), we obtain the error dynamics as 

e°(f) = -o m e°(f) + b[ki(t)x(t) + k 2 (f)r(f)] . (14) 

It is well known that the this control architecture guarantees 
asymptotic tracking x(t) — > x°(t) as t — > 00. However, the 



transient behavior of x(t) and u(t) cannot be guaranteed. The 
reason is that high frequency oscillations are generated in the 
control signal, when the adaptation rate is increased in order 
to obtain better tracking in transient. 

Next we present M-MRAC architecture. We modify the 
reference model by driving it proportionally to the tracking 
error toward the system during the transient phase, when the 
adaptive control signal u(t) is far from the ideal one because 
of the initial errors in the control gains. Thus, we design the 
adaptive control based on the modified reference model 

Xm(t) = -a m x m (t) + b m r{t) + Ae(£), x m {Q)=x 0 , (15) 

where A > 0 is a design parameter and e(£) = x(t) — x m (t). 
The e-dynamics are derived similar to the MRAC architecture 
and have the form 

e(t) = — (a m + A )e(t) + b[u(t ) — k\x{t ) + k%r(t)] . (16) 

The control signal u(t) is still given by the equation (12) as 
in the MRAC case, but the adaptive laws now are based on 
the error signal e(£) 

h(t) = -7 x(t)e(t), fci(O) = fcio 

k 2 {t) = -7 r(t)e(t), fc 2 (0) = fc 2O , (17) 

Substituting the control signal u(t) into equation (18) we 
obtain 

e(t) = ~(a m + A)e(f) + b[ki(t)x(t) + k 2 (t)r(t)\ , (18) 

where the parameter estimation errors Aq(£) and k 2 (t) are 
defined similar to MRAC case. 

In the following analysis we will also need the control 
error that is defined as u(t) = u(t) — u°(t). From the above 
constructions it follows that 

u(t) = ki(t)x(t) + k 2 (t)r(t) + k\e°{t) . (19) 

Since the ideal control signal is the best achievable signal, we 
are interested in minimizing the control error «(£), as well as 
the tracking error e°(£), both in transient and steady state by 
selecting proper values for the parameters 7 and A. This is the 
main objective of the analysis in the following sections. 

Remark 3.1: When A = 0 the proposed control architecture 
is identical with the conventional MRAC design. 

IV. Asymptotic Properties of M-MRAC 

The following theorem summarizes the asymptotic proper- 
ties of the M-MRAC architecture. 

Theorem 4.1: Let the system (8) be controlled by the M- 
MRAC scheme given by (12), (15) and (17). Then 

1) all closed-loop signals are bounded 

2) e(£) — ■> 0, e°(£) — > 0 and u(t) — > 0 as t — > 00. 

Proof: Consider the following candidate Lyapunov func- 
tion 

V{t) = ^e 2 (f) + 2- k\(t) + kl{t) . (20) 


Computing its derivative along the trajectories of the systems 
(18) and (17) we readily obtain 

V(t) = -(a m + A)e 2 (£), (21) 

which implies that e(t), ki{t), k 2 (t) £ £00 and V(t) £ Coo- 
Since the reference model (15) can be viewed as an expo- 
nentially stable system with bounded inputs r(t) and e(t), 
it follows that x m (t) £ Coo as well. Hence x(t) £ Coo- 
The boundedness of the tracking errore°(£) follows from 
the boundedness of x(t) and x°(t). Since the inclusions 
fci(f), k 2 {t) £ Coo imply that k 2 {t) £ Coo, it follows 

that u(t) £ Coo- Then u(t) £ Coo, since u°(t) £ Coo by the 
definition in (10). This completes the first part. 

Integrating the equation (21) we obtain 

(am + A) [ e 2 (r)dr = 1/(0) - V (t) . (22) 

Jo 

Since V(t) £ Coo, it follows that e(t) £ C 2 . Also, from the 
error dynamics (18) it follows that e(t) £ Coo, therefore e(t) 
is uniformly continuous. Application of Barbalat’s lemma ( 
[10], p.19) results in e(t) — > 0 as £ — » oo. Subtracting the 
equation (14) from the equation (18) we can write 

^[e(£) - e°(£)] = -a m [e(£) - e°(£)] - A e(£) , (23) 

Since e(£) — >■ 0, it follows from the equation (23) that 
[e(£) — e° (£)] — > 0, and hence e°(£) — > 0 as £ — » oo. Further, 
from the dynamic equation (8) it follows that x(t) £ Coo, 
since x(t) £ Coo and u(t) £ Coo- Therefore, x(t) is uniformly 
continuous. By assumption r(£) has a bounded derivative, 
hence is uniformly continuous. From the adaptive laws (17) 
it follows that ki(t), k 2 (t) £ Coo, therefore ki(t), k 2 (t) 
are uniformly continuous as well. Then the right hand side 
of the error dynamics (18) is uniformly continuous, therefore 
e(£) is uniformly continuous. Since e(£) has a finite limit, it 
follows from Barabalat’s lemma [5] that e(£) — > 0 as £ — > oo. 
Therefore the relationship 

lim [ki(t)x(t) + k 2 (t)r(t)] = 0 , (24) 

t — >-oo 

holds. Then it follows from the equation (19) that u(t) — * 0 
as £ — » oo. The proof is complete. ■ 

Theorem 4.1 shows that M-MRAC architecture guarantees 
not only asymptotic tracking of the modified reference model 
(15), but also the original reference model of the conventional 
MRAC. That is the asymptotic performances of both conven- 
tional MRAC and proposed M-MRAC designs are equivalent. 
However, as it is shown in the next section the transient 
performances are quite different. 

V. Transient Properties of M-MRAC 

In this section we show that the proper selection of the de- 
sign parameters 7 and A results in controllable and quantifiable 
transient performance for both output and input signals of the 
M-MRAC architecture. 



Output tracking 


To this end, we recall that Theorem 4.1 assures that the 
Lyapunov function V ( t ) is non-increasing. Therefore the fol- 
lowing chain of relationships holds 


e 2 (f) < 2V(t) < 2V (0) = - k\{ 0) + k%{ 0) 


7 L 


(25) 


which implies that 


s( 0 l < 


V7 


(26) 


where the constant a is defined as <r = y b k 2 { 0 ) + k 2 ( 0 ) 
Since (26) holds uniformly in t we conclude that 

<7 


< 


V'y' 


(27) 


Similarly, from the equation (22) the following is obtained 

a 


\e{t)\\c 2 < 


(28) 


sj (dm + A)7 

To obtain the bounds on the tracking error e°(f) we write 


the equation (23) in the following operator form 

A 


e°(s) = e(s) + 


;(s). 


(29) 


The impulse response of the transfer function G(s) = . * 

has an C\ norm of — . Therefore, it follows from (27) and 
( 1 ) that 


l|e°Wlk- < l|e(*)l U- + — l|e(t)|U,» < ( ° m + • (30) 

dm &m\ji 

On the other hand, since e(t) £ £ 2 , it follows from the 
equation (23) that e(t) — e°(f) £ £2 , hence e°(f) € £ 2 . Then, 
applying ( 2 ) and (28) we obtain 


l|e°«)lk<^J2. (31) 

a m ^n 

From the derived norm bounds in (27), (28), (30) and 
(31) it follows that the error signals e(t) and e°(t) can 
be arbitrarily decreased by increasing the adaptation rate 7 . 
However, increasing A while improves the £2 performance 
and leaves intact the Coo performance of the error signal c(t), 
also increases the norm bounds on the error signal e°(t). That 
is A cannot be selected arbitrarily and should be related to 
7 in such a way that (30) and (31) are not increased when 
increasing A. This can be achieved by setting 


a m + A = 0^7 , (32) 

where c > 0 is to be selected form the control error perspec- 
tives. In this case, we have the following norm bounds 

l|e°(t)IU« < ||e°(i)|U a < ^== , (33) 

U m CLfji -f- A 

which imply that the oscillations in the error e° (t ) are reduced 
by increasing A and 7 simultaneously according to equation 
(32), whereas the overshoot is specified by the choice of c. 

Next we analyze the transient performance of the control 
error signal u(t) defined by the equation (19). We notice 



Fig. 1. Response of MRAC to step input with 7 = 60. 


Output tracking 



Fig. 2. Response of M-MRAC to step input with 7 = 60. 


that u(t) does not explicitly depend on design parameters A 
and 7 . Instead, u(t) depend on 7 through the adaptive laws, 
and u(t) depend on A through the tracking error dynamics. 
Differentiating u(t) twice we obtain the following second 
order differential equation 

u(t) + (a m +)Xu(t) + u} 2 (t)u(t) = (34) 

7 bkip(t)e°(t ) - 7 /)(f)e(t) + n(t) + Xr 1 (t ) , 

where we denote p(t) = x 2 (t) +r 2 (t), u>(t) = \/ r ybp(t) and 
77(f) = ki(t)x(t) + k’ 2 (f)f(f) + A;*e°(f). From the results 
of the previous section it follows that all functions involved 
in the equation (34) are bounded. In particular, there exist 
positive constants 07, 0 . 2 , a 3 such that ||p(f )||£ 00 < 07, 
||p(f)IUoo < a 2 and ||r 7 (i)||£ oo < a 3 . Since the terms in 
equation (34) are also continuous functions in time, it can be 
considered as a linear equation with time varying coefficients. 
Although the equation (34) is non-autonomous, it can be still 
inferred that the adaptation rate 7 determines the frequency 
of u(t) and hence the frequency of the control signal u(t), 
since the ideal control u°(t) is in the low frequency range. 


Therefore, increasing 7 increases the oscillations in u(t) as 
it is the case for the conventional MRAC design. On the 
other hand a m + A determines the damping ratio. Therefore 
increasing A suppresses the oscillations in u(t) and hence in 
the control signal u(t). That is, by selecting a proper value for 
A the desired performance can be achieved. This is the main 
difference from the MRAC design, which results when A = 0, 
since a rn is fixed from the performance perspectives. 


Output tracking 
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Fig. 3. Response of MRAC to step input with 7 = 600. 
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Fig. 4. Response of M-MRAC to step input with 7 = 600. 


To derive the norm bound for u(t) one needs to compute 
the state transition matrix of the non-autonomous system (34). 
However, the analytic computation, even if possible, requires 
the knowledge of entire function u>(t). Since our goal is to 
derive a relationship between A and 7 , that guarantees a sub- 
optimal design, we compute an approximate bound replacing 
the function u>{t) with a constant uq. This is motivated by 
the fact that the transient behavior of u(t) mainly depends on 
the initial conditions, and by our interest in the dependence of 
the bound of u(t) on large values of the parameters A and 7 . 
That is, when the influence of the parameters A and 7 is much 
greater than the influence of other factors. Hence the constant 


b can be also ignored. Therefore, we set wq = y/'ypo- where 


Po = 


|z 0 | 
l*o I 


l£oo> 


if r(0) ± 0 
if r(0) = 0, r(t) ^ 0 


With constant ujq the equation (34) can be written in the 
following operator form 

u(s) = Gi(s)[ 7 a m p(f)e(s) + 7 fcip(f)e°(f) - 7 p(f)e(f) 
+ £t(0) — r’i(O)] + G 2 (s)[ri(f) + Tt(0)] , (35) 


where the transfer functions Gi(s) and G 2 (s) are defined in 
(4) with u> = u>o and ( = a ^ A . Applying the inequality (1) 
to each term on the right hand side of the equation (35) we 
have 


IlfiWIUo. < l|Pi(i)IU 1 [ 7 «i|^|||e°(f)|U oo + |4(0)| 

+ 7a2||e(f)|| £oo + |n(0)|] 

+ 11^2 7)|Ui [«3 + |w(0)|] ■ (36) 

Next we select 


A = 2Cov/7Po ^ a m , 


(37) 


where £0 is given in (7), and evaluate the right hand side of 
(37) term by term. Using (27), (33), (3) and (7) we obtain 


ll5i(Oll£ 1 [|«(0)| + |ri(0)|]<^[|ti(0)| + |r 1 (0)|] 
7ai|fci-H| gl (f)|U 1 l|e°(f)|U TO < 4ci r:'^ 
7a 2 |Mf)IU 1 ||e(i)|U ao < ^77 
IMOlUJas + |«(0)|] < -j=[a 3 + |u(0)|] . (38) 


Combining the terms we end up with the following approxi- 
mate bound for the u(t) 


K*)ll 


£00 7 


+ 


^ 4ciai|fci|Co_ 

<. — (7 


U rn \J P0 

2C2(«3 


2c\a2& 

PoVl 


(39) 


Vwy 


«(0)|) + 2c_i[|fi(0)| 


|n(0)|] 


Po7 


We notice that all the terms on the right hand side of (39) 
can be arbitrarily decreased by increasing the parameters 7 
except for the first term, which is independent of 7 and is 
determined by the initial estimation errors expressed by er. 

The results can be finalized in the form of the following 
theorem. 

Theorem 5.1: Let the system ( 8 ) be controlled by the M- 
MRAC scheme given by (12), (15) and (17). Then for large 
7 and A satisfying the equation (37) the following asymptotic 
bounds hold 

\\e(t)\\ Coo < o (7“ = ) , ||e(f)|U 2 < o (7“*) 

l|e°(f)|U 00 < 0(0), \\e°(t)\\c 2 <0(7-2) 

llwWIUoc < o(o-) + 0 ( 7 - 5 ) . 


(40) 


VI. Simulation Results 

For the simulation we consider the system ( 8 ) with a = 3 
and b = 2. The parameters of the reference model are set 
a m = 2 and b rn = 2. The external command is chosen to be 

1) a step function at time t = 3sec with the magnitude tq, 

2) a sinusoid of frequency 2 rad/ sec and of amplitude Tp. A 
is set to ‘2(o\/'y(xo + Tq) — a m according to equations (35) 
and (37). The response to the step command of MRAC and 
M-MRAC with 7 = 60, ro = 1 and xo = 0 are displayed in 
Figures 1 and 2 respectively. The improvement is obvious even 
with a small adaptation rate. No oscillations are observed in 
the M-MRAC response contrary to the conventional MRAC 
response. Next we increase the adaptation rate to 7 = 600. 
The corresponding performances of MRAC and M-MRAC 
are presented in Figures 3 and 4 respectively. It can be seen 
that the output tracking performance of MRAC gets better 
in the expense of high frequency control signal, while the 
output and input tracking performances of M-MRAC are both 
perfect. Further increase of 7 for the M-MRAC design results 
in the performance indistinguishable from that of the reference 
model both for input and output signals. 



Control signal 



Fig. 5. Response of M-MRAC to step input with 7 = 6000, r = 10 and 
XO = 4. 

To show the uniform performance of the M-MRAC design 
we run simulations with r o = 10 and Xq = 4, the result of 
which is displayed on Figure 5. Obviously the quality of the 
output and input tracking did not change. 

The final simulation was performed with the sinusoidal in- 
put with 7 = 6000, ro = 1 and xq = 0. The systems response 
along with the M-MRAC control signal is displayed on Figure 
5. Ones again perfect tracking is archived in both input and 
output signals. The noticeable discrepancies disappear with the 
increase of 7 . 

VII. Conclusions 

We have presented the transient and asymptotic performance 
analysis of modified reference model MRAC (M-MRAC) 
design. We showed that the desired level of tracking in 
both systems input and output signals can be achieved with 


fast adaptation provided that the design parameters are set 
according to the derived guideline. Although presented for the 
conventional MRAC design, the proposed method is suitable 
for any adaptive design technique since we are only changing 
the reference model, keeping the controller structure and the 
adaptive laws intact. 




Fig. 6. Response of M-MRAC to sinusoidal input with 7 = 6000. 
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